Role of Single Qubit Decoherence Time in Adiabatic Quantum Computation 
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We have studied numerically the evolution of an adiabatic quantum computer in the presence 
of a Markovian ohmic environment by considering Ising spin glass systems with up to 20 qubits 
independently coupled to this environment via two conjugate degrees of freedom. The required 
computation time is demonstrated to be of the same order as that for an isolated system and is not 
limited by the single-qubit decoherence time T 2 * 5 even when the minimum gap is much smaller than 
the temperature and decoherence-induced level broadening. For small minimum gap, the system 
can be described by an effective two-state model coupled only longitudinally to environment. 



Adiabatic quantum computation [l[ (AQC) is an at- 
tractive model of quantum computation (QC). It elimi- 
nates the need for precise timing of the qubit transfor- 
mations required in the gate-model computation scheme, 
and also is expected to possess some degree of fault tol- 
erance afforded by the energy gap separating the ground 
from excited states of the qubit Hamiltonian. AQC ap- 
proach is particularly appealing in the context of super- 
conducting qubits which in principle have the required 
flexibility for implementation of complicated interactions. 
In the AQC, a system starts from a readily accessible 
ground state of some initial Hamiltonian Hi and slowly 
evolves into the ground state of the final Hamiltonian Hf 
which encodes solution to the problem of interest: 



H s (t) = [l-s(t)]Hi + s(t)H f , 



(1) 



where s(t) £ [0, 1] is a monotonic function of time t. Here, 
we only consider a linear time sweep s(t) = t/tf, where tf 
is the total evolution time. Transitions out of the ground 
state can be caused by the Landau-Zener processes |2j 
at the anticrossing (s=s*), where the gap g between the 
ground state |0) and first excited state |1) goes through 
a minimum: g rn = g(s*). The probability of being in 
the ground state at the end of the adiabatic evolution is 
approximately (h — ks = 1) 



-tf/U 



to 



^9n 



as 



(2) 



To ensure large Po/, one needs tf > t a . The computation 
time is hence determined by t a and thus by g m . 

In the gate-model QC, there is no direct correspon- 
dence between the wavefunction and the instantaneous 
system Hamiltonian. The Hamiltonian is only applied 
at the time of gate operations and usually involves only 
a few qubits. The wavefunction, therefore, is strongly 
affected by the environment and is irreversibly altered 
after the decoherence time, which is typically smaller 
than the single-qubit dcphasing time T 2 *. This means 
that T 2 * imposes an upper limit on the total computa- 
tion time, unless some quantum error correction scheme 
(which requires significant resources) is utilized. This is 
not true for AQC, as the wavefunction is always very 



close to the instantaneous ground state of the system 
Hamiltonian and is consequently more stable against the 
decoherence. Qualitatively, one expects decoherence to 
drive the system's reduced density matrix towards being 
diagonal in the energy basis, which is not harmful for 
AQC but is detrimental for the gate-model QC. Such 
robustness has been demonstrated in previous studies 
S II i, H, 0, S, H G3- However, those studies have 
either used a two-state model to describe the behavior 
of a multi-level system at the anticrossing, or assumed 
noise models that are not motivated by physical imple- 
mentations. In this paper, we numerically study quan- 
tum evolution of a multi-qubit system directly without 
the two-state approximation, assuming a quite general 
and realistic coupling to environment. 

We consider a very general Hamiltonian H(t) = 
Hs(t) + Hb + Pint, which includes system, bath, and 
interaction between them, respectively. The dynam- 
ics of the total (system plus environment) density ma- 
trix is governed by the Liouville equation [lll |: p(t) = 
—i[H(t),p(t)]. The reduced density matrix for the sys- 
tem is obtained by partially tracing over the environ- 
mental degrees of freedom: ps =Tr b[p\- Let \n(t)) de- 
note the instantaneous eigenstates of the system Hamilto- 
nian: Hs(t)\n(t)) = E n {t)\n(t)) . In this basis, we define 
Pnm{t) — ( n (t)\ps{t)\m(t)) . Taking the time derivative, 
we obtain (dropping explicit time dependences) 



(n\p s \m) + (h\ps\m) + {n\p s \m). 



(3) 



We begin by focusing on the first term in ([3]) which 
is responsible for the decay processes. We treat it quasi- 
statically, assuming that the evolution of the Hamiltonian 
is much slower than the environmentally-induced decay 
rates, so that the eigenstates can be taken as time inde- 
pendent. We also assume that the effect of the system on 
the environment is so small that the bath maintains its 
equilibrium distribution at all times. Moreover, the 
bath is taken to have correlation time tb shorter than 
all the decay times of the system so that one can ap- 
ply Markovian approximation. Then, using the standard 
Bloch-Redfield formalism, one can show that [ll|, Ojl 



(n\ps\m) = -iuj nm p nm + e~ 
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(Hi,i m (0)Hi, nk (t)), 



(4) 



(5) 



(6) 



H Itnm (t) = (n\e zHBt H int (t)e 



-iH B t\ 



Here, (...) = Ttb[pb---], and summation over repeated 
indices is implied in ([5]). 

Substituting (0J into ([3]), we obtain 

Pnm = —i^nmPnm ~ E (Rnmkl ~ M nm kl) Pkl, (7) 
k,l 

where M nmkl = S nk (l\m) +S m i(n\k). The tensors M nmk i 
and Rnmkl are responsible for non-adiabatic and ther- 
mal transitions, respectively. For a time-independent 
Hamiltonian, M nmk i = 0, and becomes the Bloch- 
Redfield equations [Tl|, [l2|. The derivatives like \n) can 
be calculated numerically. It is important to ensure 
that the equation stays trace-preserving, which requires 
Re^ nm (n|m) = 0. This condition is exactly satisfied 
(even with the truncation discussed below), if we write 
(n(t)\rh(t)) = ^{(n(t+St)\ + (n(t-St)\}{\m(t+St)) - 
\m(t-St))}. 

To introduce coupling to environment, we consider a 
quite general interaction Hamiltonian 



n 

Hint -EfflM+W), 



(8) 



where aijz are the Pauli matrices of the i-th qubit, and 
Qa? are the heat-bath operators. Using ([6]), and assum- 
ing uncorrelated heat baths, we find 



•>(-) 



~E 

2 ^ 



S^(-^ k )a^ m a^ r 



fe' 



(9) 



where ^(w) = dt e^Q^ (0)) are the bath 

spectral densities and cr^ m = (Z|er£ 2 '|m). Here, we have 

neglected the imaginary parts of r>^> , as they only pro- 
duce small shifts of energies which in principle can be 
accounted for by proper renormalization. 

To model the spectral densities, we assume ohmic 



bosonic heat baths in thermal equilibrium [13| : Sa (w) 



rja we M/ W <=/(1 — e W / T ). The dimensionless coefficients 
?7a describe the strength of coupling between the qubits 




FIG. 1: (Color online). Energy spectrum for a 20 qubit in- 



stance with a very small minimum gap: g, 
Only the first 50 energy levels of the total 



,/E w 5 x 10" 
10 6 are shown. 

The dashed line in the bottom represents the ground state en- 
tanglement calculated using Meyer and Wallach entanglement 
measure (|12p . 



and environment, and u c is a cutoff frequency which we 
assume to be larger than all relevant energy scales. The 
Markovian approximation is valid as long as tb ~ l/w c 
is shorter than all the decay times and the characteristic 
time variation of the Hamiltonian. 

We now use the above model to study the evolution 
of a multi-qubit Ising system with the initial and final 
Hamiltonians given by 

§ = 45>^, do) 

i 



where Aj, hi, and Jij are dimensionless parameters and 
E is an energy scale. We consider square lattice con- 
figurations with nearest and next-nearest neighbor cou- 
plings between the qubits. The choice of only 2-qubit 
short-range interactions is motivated by the feasibility 
of experimental implementation. We generate spin glass 
instances involving 6, 9, 12, 16, and 20 qubits by ran- 
domly choosing hi and Jy from {—1,0,1} and identify- 
ing small gap instances with non-degenerate final ground 
state (see, e.g., Fig. [1]). Such instances are very rare and 
represent difficult problems; a degenerate ground state 
(multiple solutions) ensures higher probability of finding 
one of the solutions. We also choose A* = 1 for all i. 

Figure Q] shows the energy spectrum for a 20 qubit in- 
stance with a very small g m . The first two energy levels 
anticross near the middle of the evolution. We should 
mention that this is not a typical 20 qubit spectrum for 
the problem we consider. Indeed for an average prob- 
lem the minimum gap is much larger and instances with 
such a small gap are very rare. For such instances, the 
bottleneck of the adiabatic evolution is expected to be 
near the anticrossing. Moreover, since the gap is much 
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FIG. 2: (Color online). Probability of success Po/ as a func- 
tion of tf for the 20 qubit instance of Fig. [T] The solid lines 
are calculated with (j) = 0.2) and without (n = 0) coupling 
to the environment. Other parameters are E = 10 GHz, 
T — 25 mK. The case without coupling to the environment 
is compares with the pure Landau-Zener behavior using @ 
(see dashed blue line) . The excellent agreement confirms that 
for the coherent evolution, the two state model is sufficient 
to calculate the probability. The (black) dashed line, marked 
by TSM, is obtained using analytical formula (|17[) obtained 
using a two state model. It asymptotically shows the same be- 
havior as the numerical result (solid line marked by n = 0.2). 
The dotted line is numerical calculation with 77 = 0.2 6(T—g), 
which eliminates relaxation after the anticrossing (see the text 
for description). 



smaller than the typical energy separation of the levels, 
two state approximation is expected to be sufficient to 
describe the evolution. This can be tested by comparing 
a fully numerical simulation without two state approxi- 
mation with a two state model, which will be done at the 
end of this article. 

To study the evolution, we numerically integrate 
Eq. {2} starting from pg(0) — |0)(0|. For large num- 
ber of qubits, the computation becomes extremely time 
consuming because of the large number of matrix ele- 
ments in ps- However, since ps is written in the energy 
basis at all times, it is possible to significantly simplify 
the computation by truncating p$ to only the lowest few 
energy levels occupied in the course of evolution. To en- 
sure small error, we increase the number of levels kept in 
the calculation until the results do not change. Typically 
maximum 7-8 energy levels are sufficient for the type of 
evolution and instances we consider here. For very slow 
evolutions, even two states suffice to achieve acceptable 
accuracy. The numerical integration time can also be sig- 
nificantly reduced by refining the integration steps based 
on the gap size. 

Figure [2] shows the probability of staying in the ground 
state at the end of the evolution as a function of tf for 
the 20 qubit instance depicted in Fig. [T] For simplicity, 



we have chosen the same coupling to the environment 
for all qubits: rp a = t]- Notice that with the chosen pa- 
rameters, the minimum energy gap is much smaller than 
temperature (g m /T « 10 -2 ), hence thermal excitation at 
the anticrossing is expected. For a closed system (77 = 0) 
the numerics agree very well with Eq. ([2]). For an open 
system (77 = 0.2), however, the probability is enhanced 
(compared to the closed system) at small tf, while it is 
suppressed for large tf, asymptotically approaching its 
equilibrium value. Therefore, it is more efficient to run 
the system for a shorter time and repeat the process than 
to wait for a long time to attain significant probability. 
An important point is that the time scale for the prob- 
abilities to reach some non-vanishing value is almost the 
same (~ 1 ms) for all curves. This has been a generic 
property for all instances that we have studied regardless 
of the size of the gap. 

The thermally assisted behavior in short tf regime is 
the result of large relaxation after the anticrossing re- 
gion and is not expected to enhance the scaling of the 
computation [f|. To confirm this, we have repeated the 
numerical calculations, but now allowing transitions only 
in the thermal mixing region by choosing 77 = 0.2 6{T—g). 
This type of coupling coefficient only allows thermaliza- 
tion in a region with g < T and therefore eliminates the 
relaxation back to the ground state after the anticross- 
ing. The result (dotted line in Fig. [2]) shows no initial 
enhancement compared to the closed system, confirming 
the above statement. 

We now compare the numerically calculated compu- 
tation time with the single-qubit decoherence times. If 
the qubits are uncoupled (Jij—0), the single-qubit de- 
coherence rates, in weak coupling limit, are given by 
1/T 2 * - S&(0) ~ 77T. For the parameters of Fig. 
T 2 *~10 ns, which is typical for solid-state qubits. This 
decoherence time is five orders of magnitude smaller than 
the computation time (~1 ms) for the problem of Fig. [2] 
Therefore, unlike the gate model QC, in AQC, the com- 
putation time is not limited by the single-qubit decoher- 
ence time. 

It should be noted that the qubits will go through an 
entangled state during the evolution. To demonstrate 
that, we have displayed in Fig. [1] the ground state en- 
tanglement (dashed line) calculated using the measure 
originally proposed by Meyer and Wallach [141 ]: 

1 - 

Q(|V)) = -^2(l-Tr[ (0 2 fe ]), (12) 
fc=l 

where pk = Ttj^klip) (ip\ is obtained by partially tracing 
over all qubits except the k-th one. The ground-state en- 
tanglement becomes nonzero in the first half of the evolu- 
tion with a very sharp peak at the anticrossing. The mea- 
sure (|12p. however, does not describe quantitatively the 
actual (mixed-state) entanglement at that point, since 
it does not account for thermal mixing. Unfortunately, 
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no practical mixed-state entanglement measure exists for 
more than two qubits [l5T | . The entanglement will not be 
destroyed by the environment (except maybe at the anti- 
crossing) as long as the system dominantly populates the 
ground state. Therefore, despite the fact that the evolu- 
tion time is far beyond the qubits' dephasing time, the 
system still preserves its quantum mechanical behavior 
throughout the evolution. Once again, this is in contrast 
to what is expected in the gate-model QC. Qualitative 
demonstration of the non-vanishing ground state entan- 
glement during the evolution is also important as it shows 
that the evolution cannot be described efficiently in only 
classical terms. 

The numerical method presented here is valid only for 
a Markovian environment. Most environments, however, 
especially in superconducting systems, are not Markovian 
and have a significant amount of low frequency noise. 
In Ref. [l(J, we have used the two-state model (TSM) 
approximation to study the effect of a non-Markovian 
environment on AQC. For the rest of this paper, we focus 
on showing that the TSM is adequate for the description 
of the AQC performance in the small-gap regime. For the 
environment to be able to cause transitions out of the 
ground state, the interaction Hamiltonian should have 
nonzero matrix elements between the ground state and 
the target state. We introduce 



(13) 



which give the rms values of the matrix elements of Pauli 
matrices a4 between the lowest two states. They repre- 
sent some average behavior of the corresponding matrix 
elements and are closely related to relaxation between the 
two levels. Especially, if all the qubits have the same cou- 
pling r\ to the environment, the relaxation rate between 
the two states is given by 




7 = n(M 2 x + M 2 z )S(u w ), 



(14) 



where S(uj) is the symmetrized spectral density of the 
(uncorrelated) baths. Figure[3](a) displays M x and M z as 
a function of s for the 20-qubit system of Fig. [TJ Except 
for the initial region, they both show the same behavior: 
a sharp peak at the anticrossing, with a width propor- 
tional to g m , followed by a vanishingly small value. For 
small T, the excitation from the ground state will be sup- 
pressed everywhere except near the anticrossing, where 
g<T. 

The transitions at the anticrossing can be described by 
an effective two-state Hamiltonian: 



H 



TSM 
S 

TSM 



~(g m T x + er x ), 
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TSM 



T X + QJ SM T Z ), 



(15) 



where r Q 
space, e = 



are the Pauli matrices in the two-state sub- 
; 2E(s — s*), with s* being the position of the 



0.8 
0.6 

s 

0.4 
0.2 




(a) 
















) 0.2 


0.4 


0.6 


0.8 1 




s 








0.52 



FIG. 3: (Color online), (a) The rms value of the matrix el- 
ements of the Pauli matrices between the first two states, as 
define in (|13|l . Solid (blue) line is M z , dashed-dotted (black) 
line is M x . We have also plotted the matrix element (Mj SM ) 
of t z between the two energy levels of a two state model de- 
scribed by Hamiltonian (|15|l , (b) The same curves zoomed 
near the anticrossing show qualitative agreement (except for 
a prefactor) between all three curves. 



anticrossing and E an energy scale (in our case, close 
to E) characterizing the anticrossing. We introduce the 
matrix elements 



Mj SM 

M TSM 



|t s |1)| = | cost 
ItJDI = Isin6 



(16) 



which are responsible for the transitions between the two 
levels, where tan0 = g m /e. Except for an overall fac- 
tor, AfJ SM has the same shape as both M x , z (see the 
inset of Fig. [3]), while Mj SM is completely different; it 
has a sharp dip at e = where it vanishes. Thus, in 
order for the effective two-state model to give the same 
result as the full system, it needs to couple to the envi- 
ronment only via t z (i.e., Qj SM = 0). In other words, a 
generic single-qubit coupling to the environment reduces 
to only longitudinal coupling in the effective TSM. Be- 
sides the numerical agreement, the fact that Qj SM = 
in the two state model has a deep physical meaning. As 
is clear in (fTrj)) . A1 X SM w 1 everywhere except for a very 
small region near the anticrossing. In the presence of 
a non- vanishing Q X SM , the relaxation rate will be very 
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large almost everywhere and therefore the system should 
relax to the ground state in constant time even when 
s = 1. This is obviously unphysical and especially not 
expected for spin glasses for which the relaxation time to 
the ground state is extremely long. In the appendix of 
Ref. the effective Hamiltonian (TT5)) is systematically 
derived for the case of adiabatic Grover search problem 
[IB]. For that problem, one finds Q^sm = o(\/y/N), but 
qTSM _ 0(1), and therefore in large N limit the former 
vanishes in agreement with our physical expectation. No- 
tice that in this case, the relaxation rate due to coupling 
of the bath to t x is 7 cx (Qj SM ) 2 = 0(l/N), therefore 
solving the problem merely based on relaxation leads to 
a computation time tf = O(N), which is the complexity 
of classical computation. 

From the TSM Hamiltonian ([15]) with this type of lon- 
gitudinal coupling, the success probability in the large-T 
limit is [lUcl 



pTSM 
0/ 



-2t f /t 



)• 



(17) 



This formula is also plotted in Fig. [2l The qualitative 
asymptotic agreement with other numerical curves in the 
figure indicates that most of the transition occur in the 
small-gap region (g <C T) where the TSM is adequate for 
their description. 

To summarize, by studying spin glass instances of up 
to 20 qubits (only one is illustrated), we have explicitly 
demonstrated that the computation time in AQC can be 
much longer than the single-qubit decoherence time T 2 *. 
In the case of small minimum gap (i.e., hard instances), 
effective two-state model with only the longitudinal cou- 
pling to environment describes transitions at the anti- 
crossing. The numerical results also show that the com- 
putation time scale is unaffected by ohmic environment. 
This conclusion can not be understood directly as sup- 
pression of transitions to the excited states by the energy 
gap, since in the chosen instances it was much smaller at 
the anticrossing than the temperature and decoherence 
strength. Rather, it arises from the balance of transitions 



between the two lowest levels. It should be emphasized 
that these results were obtained under the assumption of 
weak coupling to the environment, for which the discrete 
energy structure of H$ is mostly preserved. In the case of 
strong coupling, the interaction Hamiltonian would dom- 
inate, and both the method and the conclusions of this 
work would not hold. 
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